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Abstract 

From recent CMB and Large Scale Structure observations the value of the equation 
of state of dark energy, assuming it to be constant in time, is constrained to be 
— 1.3 < uu < —0.7 at the 95% confidence level: consistent with dark energy being a 
classical vacuum term. Here we describe two novel and independent methods, sen- 
sitive to different systematics, that give the same value for w and similar confidence 
regions. This suggests that systematics are not an issue in current determinations of 
w. The first method yields a measurement of w that relies on the minimum number 
of model-dependent parameters; the second method is a non-parametric measure- 
ment of the time dependence of w(z). We also present a method to statistically 
determine the edge of a distribution. 
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1 Introduction 



It is n ow established that the universe is ac c elerating (e.g.. iDunlop et al. 
f)l996MSpinrad et all |l997h : iRiess et all (|l998MPerlmntter et all (|l999h : the 
next logical step is to unveil the nature of the accelerating force, for example, 
by determining its equation of state {p = w(z)p), where w(z) need not to be 
constant and can vary with redshift. 



During the past year we have seen a significant improvement in the accu- 
racy of th e equation of state, w measurements ( assum i ng w to be con s tant i n 
time, e.g. ISpergel et all (|2003h : ICaldwell et all (|2003lb I.Timenez et all (|2003h : 
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Tonrv et alJ ((20031 ) and also Verde and Melchiorri articles in these proceed- 
ings). All these methods are sensitive to different systematics, yet the results 
are re markably consistent: —1.3 < w < —0.7 at 95% confidence (ppergel et al. 



2003). While this measurement will soon become even tighter with 2yr WMAP 
data and SDSS galaxy and Ly— a forest power spectrum release, it still assumes 
w to be constant in time. A significant challenge will be to accurately mea- 
sure w(z). There h as been a number of p a pers di scussing different method s 



is been a number ot p a pers discussin g dirterent metnod s 
to mea sure w ( e.g. I Alcock and Paczvnskil (I1979h: Turner and White (Il997h: 
FeTall (ll998h:lGarnavich et al l (Il998l ) : iBirkinshawl (Il999l): lEfstathiou 



Caldwc 



1999); 



Hui (1999): Newman and Davi 



Haiman et al 



200l[). 



Huterer and Turner 



Wang and Garnavichl (12001): Weller anc 


Albrec 


ht|( 


2001);] 


3accigaluDi et al.l 


Hutererl(|2002h:lHul( 2002b!) : Kuiat et al. 


(20021): 


Hu 


d2002a 


) : Lima and Alcanizl 



(2001); 



(2002); 



()2002f l) while less attention has been devoted to the more challenging task of 
measuring w(z) without parameterizations. 



Below, we describe two novel methods: a) a measurement of w that relies 
on the minimum number of model dependent parameters. This method uses 
the determination of the position of the first acoustic peak in the CMB an- 
gular power spectrum (£i) and the absolute ages of galactic globular clusters 
(GCs). b) a non-parametric measurement of the time dependence of w(z). This 
method is based on th e relative ages of st e llar p o pulations. These tech niques 



metnod is based on tn e relative ages oi st e llar p o pulations, inese teem 
have been described in Jimenez and Loebl ( 2002^ 1: Ijimenez et al. (2003) 



2 w from the t\— GCs ages method 



If the universe is assumed to be flat, the position of the first acoustic peak {l\ 
in the standard spherical harmonics notation) depends primarily on the age of 
the universe and on the effective value 1 of w (jCaldwell et al. Ul998l Eu" et all 
200lt iKnox et all 1200 it ICaldwell et "all . 120031) . As noted by these authors, for 
a fixed w value, a change in the physical density parameter ui m = Q m h 2 that 
keeps the characteristic angular scale of the first acoustic peak fixed will also 
leave the age approximately unchanged. Thus an independent estimate of the 
absolute age of the universe at z = combined with a measurement of l\ 
yields an estimate of w largely independent of other cosmological parameters 
(e.g., Q m and h). 



This can be better understood by considering that for a constant w(z), the 



1 I.e., the average value over redshift. We use w(z) to indicate when w is allowed 
to vary over time. 
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age of a flat universe is given by 



t = h, 1 I (i + z)- 1 \n xad (i + z f + n m (i + zf + o A (i + z) 3 ^ 1+w) 



-1/2 



dz .(1) 



The position of the first acoustic peak is fixed by the quantity 

9a = r s {a dec )/d A {a dec ) , (2) 



where a dec is the scale factor at decoupling, r s (ad ec ) is the sound horizon at 
decoupling, and d,A(a>dec) is the angular diameter distance at decoupling. For 
a flat universe, 



a de 



r s (a dec ) 



1 + ^) (Hax 1 - 3 '" + fi m x + fi rad ) 



-1/2 



cfx (3) 



and 



dAi^de 



c 



Max 1 - 3 ™ + fi m x + Q 



rati 



-1/2 



(4) 



See Verde et al. (2003) for more details. We use equation s 1 to 4 and fix 9a 



to be consistent with the WMAP determination ~ 220; Page et al. . 20031 ) 



Figure 1 shows the allowed region in the age-u> plane obtained for VL^h 2 = 0.02, 
0.05 < Q m < 0.4, and 0.5 < h < 0.9. It is clear from the plot that an 
independent and accurate age determination of the universe will provide a 
measurement of w. Thus by using only the WMAP observation of t\ and an 
independent estimate of the age of the universe, one can place constraints 
on w, largely independent of other cosmological parameters. Note that this 
method is sensitive to a different redshift weighting than supernova and CMB 
measurements. 



The ages of the oldest globular clusters (GCs) provide a lower limit to the 
total age of the universe. Since nume rous star-forming gal axies have now been 
observed up to redshift z = 6.6 (e.g., Kodaira et al. . 20031 ) and the oldest GCs 
contain the oldest stellar populations in galaxi es, it is reasonable to a ssume 
that GCs too have formed at redshift z > 6. Krauss and Chabover (2003) 
perform the most careful analysis to date of the effects of systematics in GC 
age determinations. They estimate the age of the oldest Galactic GCs using 
the main-sequence turn-off luminosity and evaluate the errors with Monte- 
Carlo techniques, paying careful attention to uncertainties in the distance, 
to systematics, and to model uncertainties. They find a best-fit age for the 
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Fig. 1. The region of parameter space allowed if we fix 9a (which in turns fixes the 
position of the first CMB acoustic peak) and uji = Q,bh 2 = 0.02, weakly constraining 
0.05 < fl m < 0.45 and 0.5 < h < 0.9. 

oldest Galactic GCs of 12.511:2 Gyr (95% confidence limits). We find that their 
probability distribution for the oldest GC age, can be accurately described by 



P(t) 



A 

<t - T) 



cxp 



In [(t - T)/m] 
2o~ 2 



(5) 



where t denotes the age of the oldest GCs, A = 1.466, a = 0.25, T = 6.5 
and m = 5.9. Since the oldest GCs likely formed at z > 5, for all reasonable 
cosmologies we only need add about 0.3 Gyr to the GC ages to obtain an 
estimate of the age of the universe. Only the oldest (low metallicity) GCs 
should be used i n the a ge estimate, for which there is no age spread (see e.g. 
IRosenberg et~aT1 (|l99flj) . which finds no age spread for GCs with metallicities 
[Fe/H] < -1.2). 

In order to constrain w, we can now assume a flat universe and Monte Carlo 
simulate the distribution of i\ subject to only weak priors on the ot her cosmo- 
logical parameters. For different values of w. I Jimenez et al. (2003) generated 



10 5 models, randomly sampling the cosmological parameters Q m , flbh 2 and h 
with uniform priors, 0.5 < h< 0.9, 0.05 < ft m < 0.45, and 0.01 < Vt h h 2 < 0.03. 
The next step is to impose an age of the universe constraint by randomly sam- 
pling these models so that the age of the universe has a p robability distribution 
whose shape is given by equation 5, but offset by 0.3 Gyr. Jimenez et al. (l2003| ) 



used the publicly available code CMBFAST ( Seliak and ZaldarriagaL 1996i l to 
compute the resulting distribution of t\. This is shown in Figure 2 where the 
two solid lines are the 68% and 90% joint confidence levels. 

As expected, the age alone does not constrain w; there is a degeneracy between 
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Fig. 2. Confidence regions (68 and 90%) for the value of w are shown by the filled 
grey areas (marginalized over the position of the first acoustic peak in the CMB, 
The solid lines denote contours which show the constraints in w using the derived 
age of the oldest Galactic globular clusters only. The dash and dotted lines are the 
the position of the first acoustic peak and confidence region as measured by WMAP. 
Note that the position of the first peak greatly constrains the value of w. 

w and l\. If we now use the observed position of the first acou stic peak as 



recen tly measured, in a model independent way, from WMAP (jPage et al 
2003), we can constrain w with high accuracy. The filled contours in Figure 2 
are marginalized values for w at the 68% and 90% confidence levels. Thus we 
find w < -0.8 (w < -0.67) at the 68% (90%) confidence level. 



This determination depends solely on the GC determination of the age of the 
universe and on the observed position of the first acoustic peak in the CMB 
power spectrum. Th is constraint is slightly less stringent than that obtained by 
Spergel et al.l (j^003) from a joint likelihood analysis of WMAP with six exter- 
nal data sets ( WMAP + CBI + ACBAR + 2dFGRS + Lymana forest power 
spectrum + Type la supernovae + H constraint from the HST key project), 
but is tighter than the CMB-only ( WMAP + CBI + ACBAR) determination 
and comparable to the WMAP + ACBAR + CBI + HST constraint. 



3 Measuring a variable equation of state for the dark energy 



The popular approach for measuring w(z) uses its effect on the luminosity 
distance of sources. In particular, the proposal for the Supernova/Acceleration 
Probe (SNAP) mission 2 suggests to monitor ~ 2000 Type la supernovae 

2 http://snap.lbl.gov/ 
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across the sky per year and determine their luminosity distances up to a 
redshift z ~ 1.5 with high precision. However, the sensitivity of the luminosity 
distance to the redshift evolution of w(z) is compromised by its integral nature 
(Maor et al. 2001), 

o 

r dt 

d L = (l + z) J (l + z')—dz', (6) 



where t(z) is the age of the Universe at a redshift z which depends on w(z). 



Jimenez and Loe"bl (|2002l ) proposed an alternative method that offers a much 



better sensitivity to w(z) since it measures the integrand of equation (6) di- 
rectly. Any such method must rely on a clock that dates the variation in the 
age of the Universe with redshift. The clock is provided by spectroscopic dat- 
ing of galaxy ages. Based on measurements of the age difference, At, between 
two passively-evolving galaxies that formed at the same time but are sepa- 
rated by a small redshift interval Az, one can infer the derivative, (dz/dt), 
from the ratio {Az/ At). The statistical significance of the measurement can 
be improved by selecting fair samples of passively-evolving galaxies at the 
two redshifts and by comparing the upper cut-off in their age distributions. 
All selected galaxies need to have similar metallicities and low star formation 
rates (i.e. a red color), so that the average age of their stars would far exceed 
the age difference between the two galaxy samples, At. 

This differential age method is much more reliable than a method based on an 
absolute age determination for galaxies (e.g., Dunlop et al. 1996; Alcaniz Sz 
Lima 2001; Stockton 2001). As demonstrated in the case of globular clusters, 
absolute stellar ages are more vulnerable to systematic uncertainties than 
relative ages (Stetson, Vandenberg & Bolte 1996). Moreover, absolute galaxy 
ages can only provide a lower limit to the age of the Universe and only place 
weak constraints on the possible histories of w(z). 

Consider a flat universe composed of matter and dark energy with an equa- 
tion of state pq = w(z)pq. The Hubble parameter is H 2 = Hq[pt(z) / pt(0)}. 
Here, the subscripts Q, m, and T refer to the dark energy, the matter, or 
the total sum of the two, respectively. Assuming further that the matter is 
non-relativistic (i.e. effectively pressureless) , we get 



r dz' 

On(0) + n Q (0) exp 3 J ——w{z) 



1/2 

, (7) 



where we have used the energy conservation equation for the dark energy, 
Pq = —3H(1 + w(z))pq (Maor et al. 2000). Thus, (dz/dt) is related to the 
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equation of state of the dark energy through one integration only, while the 
luminosity distance in equation (6) is given by an integral of the inverse of 
equation (7), namely through two integrations. By differentiating (dz/dt) with 
respect to t we find, 

-^ m (0) (l + z) A w{z), (8) 



rr-2 d Z 



[H \dz/dt)Y 

(T+z) 



2 + 2 W(Z) 



which depends explicitly on w(z) without any integrations. Thus, the second 
derivative of redshift with respect to cosmic time measures w(z) directly. While 
it is possible to find significantly different redshift histories of w(z) for which 
the evolution of g?l is similar, this cannot be done for (d 2 z/dt 2 ). 



To obtain dz/dt, Ijimenez and Loebl (|2002j ) proposed to use the old envelope 
of the age-redshift relation of E/SO galaxies. They assume that E/SO at dif- 
ferent redshifts are drawn from the same parent population wit h the bulk of 
their stellar populations fo rmed at relatively high redshift (e.g.. iBower et aL 



1992; IStanford et all Il998 h . Then, at relatively low redshift, they are evolving 
passively and may be used as "cosmic chronometers" . 



Recently, I Jimenez et al. (2003) applied this method using high-quality spectra 
of a sample of old stellar populations covering the largest possible range in 
redshift. To this aim they combined the following data sets: (i) the luminous 
red gal axy (LRG) sample from the Sloan Digital Sky Survey (SDSS) early data 
release feisenstein et al. . 200lh : (ii) the sample of field early-type galaxies from 
Treu et all (ll999ll200lL 12002. hereafter the Treu et al. sample); (iii) a sample of 
red galaxies in the galaxy cluster M S1054— 0321 at z = 0.83; and (iv) the two 
radio galaxies 53W 091 and 53W069 (|Dunlop et al.Ul99d Ispinrad et al!ll997 : 
Nolan et all 120081 Dey et al., in preparation). They obtained the age of the 
domina nt stellar population i n the galaxies by fitting single stellar population 
models ( Jimenez et al. . 1998| ) to the observed spectrum. 



Figure 3 shows the derived single-burst equivalent ages of the galaxies as 
a function of redshift. The circles correspond to galaxies in the SDSS LRG 
sample, triangles refer to the Treu et al. sample, diamonds are galaxies in 
MS1054-0321, and crosses are 53W091 and 53W069. For clarity, the SDSS 
LRG points have not been plotted for ages < 7 Gyr. Typical errors on the 
ages of LRG galaxies are 10% and are not plotted. An age-redshift relation 
(i.e., an "edge" or "envelope" of the galaxy distribution in the a ge-redshift 
plane) is apparent from z = to z = 1.5. Ijimenez et aD (|2003l) also show 
that small recent episodes of start formation do not significantly bias the age 
determination. 



What are the observational limits that can be imposed on w(z)l As a con- 
sistency check we note that the age at z = obtained with this method is 
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Fig. 3. Top panel: Age-redshift envelope obtained for the galaxies studied in 



Jimenez et al.l (|2003). Open circles correspond to galaxies from the SDSS LRG sam- 
ple, triangles to the Treu et al. sample, diamonds to MS1054— 0321, and crosses are 
53W091 and 53W069. A clear trend is present: galaxies age as the redshift decreases. 
The overall shape of this trend is in fair agreement with theoretical expectations 
for a ACDM (solid line) with Q, m = 0.27 and H = 71 km s" 1 Mpc" 1 . The dotted 
line represent an alternative model disfavored by the data: w(z) = —2 for z > 1 
and then grows linearly from w = —2 at z = 1 to w = at z = 0. Bottom panel: 
histogram along the age axis of top panel for the redshift range < z < 0.04 and 
0.08 < z < 0.12. The clear shift between the two histograms is a measurement of 
dz/dt and therefore allows us to measure Hq. 



in good agreement with GC ages and that the Treu et al. galaxy ages agree 
with those from of SDSS LRG sample where they overlap. The solid line in 
Figure 3 corresponds to the age-redshift relation for a flat, reference ACDM 
model: tt m = 0.27, h = 71, and w = -1 (WMAP best fit model). This is 
consistent with the observed age-redshift envelope, and it seems to indicate 
that we live in a universe with a classical vacuum energy density and that, 
on average, stars in the oldest galaxies formed about 0.7 Gyr after the Big 
Bang (e.g., at redshifts z ~ 7 — 10). To make the observed age-redshift rela- 
tion consistent with a model with widely different w(z) behavior, one would 
have to infer that early-type galaxies at different redshifts have very different 
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formation epochs for their stellar populations. For example, the dotted line in 
Figure 3 indicates the age-redshift relation for a passively evolving population 
in a universe where w(z) = —2 for z > 1 and grows linearly from w = —2 at 
z = 1 to w = at z = 0. For this model to work, we would have to conclude 
that the oldest galaxies at z ~ 1.5 formed about 0.7 Gyr after the Big Bang, 
but that the oldest galaxies at z ~ formed ~ 3 Gyr after the Big Bang. 
Since the oldest galaxies in the LRG sample are as old as the oldest GCs, we 
consider this to be an unlikely explanation. Furthermore, not only would this 
scenario require that all local galaxies either formed their stars more recently 
than galaxies at high redshift or had recent episodes of star formation, but 
would also require a remarkable fine-tuning for them all to end up with ex- 
actly the same age. We therefore conclude that this extreme model for w(z) 
is unlikely given the data. 



Unfortunately, the small number of galaxies in ljimenez et al. (2003) samples 
at z > 0.2 did not allow them to compute dz/dt with enough accuracy to 
constrain w(z), and therefore this measurement will have to await better data. 
However, the z < 0.2 region is well populated by LRG galaxies (Figure 3) and 
they determined dt/dz at z ~ 0. This is illustrated in the bottom panel of 
Figure 3 by the clear shift in the upper envelope of the age histogram between 
the two redshift ranges < z < 0.04 and 0.08 < z < 0.12. This can be used 
to determine Hq, which relies upon determining the "edge" of the galaxy age 
distribution in different redshift intervals. 



If the ages of the galaxies were known with infinite accuracy, for each galaxy j 
at redshift Zi, one could associate an age a;. The probability for the age of the 
oldest stellar population in galaxy i; P(t|oj), would be given by a step function 
which jumps from to 1 at <2j. An age-redshift relation edge could then be 
obtained by dividing the galaxy sample into suitably-large redshift bins and 
multiplying the P(t\a,i) for all the galaxies in each bin. 



In practice, the age of each galaxy is measured with some error <faj. We thus 
assume P(t\a,i) = 1 if t > a, + S n , and InP fa,) = — x 2 where x = (a, + 5 ai — 
t)/(y/25 ai ) otherwise. Jimenez et al. ( 2003| ) divided the z < 0.2 portion of 



the LRG sample into 51 redshift bins. For each bin, they obtained P(z,t) by 
multiplying the P(t\a,i) of the galaxies in that bin. They define the "edge" of 
the distribution t(z) to be where In P(z, t) drops by 0.5 from its maximum (i.e. 
A InP = 0.5 and associate an error to this determination 8h z \ given by t 2 (z) — 
t(z) where £2(2) corresponds to where lnP(^, t) drops by 2 from its maximum 
(i.e. AlnP = 2, this approximately corresponds to the 68% confidence level). 
This procedure makes the determination of the "edge" less sensitive to the 
outliers (see also Ravchaudhurv et al. ( 1997i l). Figure 4 presents the resulting 
t(z) relation. 
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Fig. 4. The binned, low-redshift age-redshift relation, derived from the SDSS LRG 
sample. The solid line is a best-fit to the edge. 
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Fig. 5. The value for Hq derived from Figure 4 as a function of redshift cut-off for 
computing the slope. Note that it is very weakly dependent on the value of this 
cut-off. The solid circle is the adopted value. The shaded are corresponds to the la 
confidence region from the Hubble Key Project. For z > 0.18 the probability of the 
fit being a straight line is P> x z < 0.1, while at z = 0.17 is 0.32. 




For z < z max , I.Timenez et al.l (j2003[ ) fit t(z) with a straight line whose slope 



dt/dz is related to the Hubble constant at an effective redshift by H(z e s) 
— 1+ * - This fit is performed by standard x 2 minimization. They also 

compute P> x 2, the probability of obtaining equal or greater value of the re- 
duced x 2 if t ne t(z) points were truly lying on a straight line. Values of 
P> x 2 < 0.1 means that a straight line is not a good fit to the points. 
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For w = -1, H = H(z cS ) [Q m (l + z cS ) 3 + fi A ] 2 where fi m = 0.27, fi A = 0.7. 
Figure 5 shows how the H measurement depends on z max . For z max > 0.17 not 
only the value for H drifts but also a straight line is not a good fit to the points 
(P> x 2 suddenly drops below 0.1). For their H determination, z max = 0.17, 
P> x 2 = 0.32, z cff = 0.09, and the correction from H(z e s) to H is a 4% effect. 
They obtained Hq = 69 ± 12 km s _1 Mpc -1 . Which is in goo d agreement with 



the H ubble Key Project measurement (h = 0.72±0.03±0.07: lFreedman et al 



2001), with the valued derived from the joint likelihood anal ysis of WMAP 



+ 2d FGRS + Lyman-a forest power spectrum (h = 0.7llo;og; ISpergel et al. 



2003 ), with gravitational lens t ime delay determinations (h = 0.59lo!o7 ^ 



0.03; Treu and Koopmans l. 20021) . and Sunyaev-Zeldovich measurements (h 



n fin+°- 04+0 - 13 - 

u - uu 0. 04-0.18 ) 



R.eese et alll2002h . 



This check at z = gives some confidence about the possibility of measuring 
w(z) when larger samples of elliptical galaxies at higher redshifts become 
available in the near future. DEEP2 and the full SDSS catalog will provide 
excellent datasets for doing this. 

I warmly thank my collaborators Avi Loeb, Dan Stern, Tommaso Treu and 
Licia Verde. 
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